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ABSTRACT 

Spectral distortions of the CMB have recently experienced an increased interest. One 
of the inevitable distortion signals of our cosmological concordance model is created 
by the cosmological recombination process, just a little before photons last scatter 
at redshift 2 ~ 1100. These cosmological recombination lines, emitted by the hydro¬ 
gen and helium plasma, should still be observable as tiny deviation from the CMB 
blackbody spectrum in the cm-dm spectral bands. In this paper, we present a forecast 
for the detectability of the recombination signal with future satellite experiments. We 
argue that serious consideration for future CMB experiments in space should be given 
to probing spectral distortions and, in particular, the recombination line signals. The 
cosmological recombination radiation not only allows determination of standard cos¬ 
mological parameters, but also provides a direct observational confirmation for one of 
the key ingredients of our cosmological model: the cosmological recombination history. 
We show that, with present technology, such experiments are futuristic but feasible. 
The potential rewards won by opening this new window to the very early universe 
could be considerable. 
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1 INTRODUCTION 


The cosmic microwave background (CMB) provides one 
of the cleanest sources of information about the Universe 
in which we live. In particular, the CMB temperature 
and polarisation anisotropies have allowed us to pin down 
the key cosmological parameters with unprecedented 
precis ion (iBennett et al.l l2003l : IPlanck Collaboration et al.l 
l2014al ). and we are presently witnessing the final stages 
in the analysis of Planck data dPlanck Collaboration et al.l 
l2015bl l. Most of today’s experimental effort is going into 
a detection of primordial B-modes, exploiting the curl 
polarisation p atterns sourced by gravity waves created dur¬ 
ing inflation (iKamionkowski. Kosowskv fc StebbinsI 199j 


Se ^k fc Zaldarriaeal 1 19971 : IKamionkowski fc Kosowskv! 

Many balloon or ground-based experiments, 
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such as SPIDER, BICEP2, K eck Array, Simons Array . 


CLASS, etc, (e.g.. ICrill et al.l 2008: lEimer fc et al.l l20ld: 


BICEP2 Collaboration et al.l 


2014 


Staniszewski et al] 


2012h . are currently observing or coming online. These 


measurements will eventually exhaust all the information 
about the primordial Universe to be gained from the CMB 
anisotropies. 

The next frontier in CM B measurements is t he de¬ 
tection of spectral distortions (|Silk fc Chlubal f2014l ~l. This 
is a field that has not changed since COBE/FIRAS set 
upper lim its on the ti and y parameters more than 20 
years ago dFixsen et al.lflO^ L In the standard cosmological 
model , chemical potential distortions dSunvaev fc Zeldovichl 
Il970ch are expected at a level of about 10 ^ of the FI- 
RAS upper limits i.e., ~ 10“®) as a consequence of 

the damping of primordial a diabatic density fluctuations 
I Su^amA_Zeldoyichj_ 1970^ Da^ 199ll : IHu. Scott fc SilkI 
Il994al : Chluba. Khatri fc Sunvaevl2012l 'l . A detection of this 
signal would probe the redshift range between thermalisa- 


© 0000 RAS 























































2 Desjacques et al. 


tion ( 2 : ~ 2 X 10®) and recombination (2 ~ 10®). While very 
interesting constraints cou ld be derived for non-standard 
inflat ion scenarios (e.g., IChluba. Erickcek fc Ben-DavanI 
l2012h . since the distortion is cause by an integrated ef¬ 
fect of energy injection in the early Universe, it would 
be challenging to distinguish a detection from more ex¬ 
otic sources of energy, such as damping of small-scale non- 
Gauss ian density or blue-til t ed primordial density fluctua- 


IGhluba & Grin 

20131. or even dark matter self-annihilation 

IGhluba 2013a 

Ghluba & Jeonell2014h. 

A truly powerful probe of spectral distortions must 
carry redshift information. While the hybrid non-/x/non- 

V distortion | 

llhluba & SunvaevI l2012l: iKhatri & SunvaevI 

I 2 OI 2 I: IGhluba 

range (10"* < 

probe, capable 
precedentedly 
helium recom 

l2013bll is limited to a narrow redshift 
2 < 3 X 10®), a far more powerful 

of probing the Universe directly at un- 
high redshift, would be the hydrogen and 
aination line snectrum dDubrovichl 1975: 

Rvbicki delP Antonio! 19931: Dubrovich Sz StolvarovI 1995. 


in the crown of any future spectral distortion experiment, 
capable of verifying the recombination history of the Uni- 
verse and measuring the primordial helium frac tion (e.g., 
IChluba &L Sunvaevl2008l : [Sunvaev fc Chluball2009l l. It would 
simultaneously provide a unique calibration template for 


u-distortion dZeldovich & Sunvaev 

I 1969 I 1 caused b-v 

^ reion- 

ization and structure formation ( 

Hu. Scott & SilkI 

Il994bl: 


mentioned dissipation signal. A cosmic recombination line 
probe would require extensive frequency coverage over GHz 
to THz frequencies and exquisite sensitivity. It would open 
a new window on the Universe 380,000 years after the Big 
Bang (ISunvaev k, Ghlu^l2009ll . 

Here, we explore the feasibility of such a probe for a 
future space experiment. We will consider the following sit¬ 
uation: low- and high-frequency channels are used to remove 
synchrotron radiation and thermal dust emission, respec¬ 
tively. At the remaining available frequencies, we can use a 
template for the recombination spectrum and cross-correlate 
it with the data. This matched-filtering approach is partic¬ 
ularly well suited for our purpose as we are mainly inter¬ 
ested in the detection level of the recombination spectrum. 
It will naturally take into account frequency correlations in 
the template. At the Fisher matrix level, this is equivalent 
to a forecast for the overall amplitude of the recombination 
spectrum template. 

The paper is organised as follows. We begin with a de¬ 
scription of our model for the sky signal in l|2j with particu¬ 
lar emphasis on the cosmic infrared background (GIB). We 
discuss experimental setups and present forecasts in M We 
conclude in [jd] We adopt a concordance, flat AGDM model 
in all illustrative calculations. 


2 MODELLING THE RESIDUAL INTENSITY 

We are interested in modelling the specihc intensity (or 
brightness), Imaiv), that remains after the subtraction of 
the average CMB blackbody spectrum as a function of fre¬ 
quency. Assuming that foregrounds such as dust etc. have 


been subtracted, this residual intensity is a sum of GMB 
spectral distortions, the GIB and noise induced by the inci¬ 
dent radiation and the detector. We will now describe each 
of these components in more detail. 


2.1 Spectral distortions of the CMB 

Upon subtracting the estimated GMB monopole, we are left 
with primordial /r- and ^/-distortions, a temperature shift, 
which arises from our imperfect knowledge of the GMB tem¬ 
perature, and the cosmological recombination radiation. We 
shall neglect the residual (non-/r, non-y) distortion signal re¬ 
lated to the precise time-de pendence of the energy release 
process (e.g., lGhluball2013bl l. which can add extra low in¬ 
tensity features to the broad primordial distortion. 

For the signals, we use the dehnitions commonly 
adopted in the treatmen t of GMB spectral distortions (e.g ., 
iGhluba fc Sunvaevll2012l : iGhluba. Khatri &: Sunvaevll2()12 1 . 

For the /r-distortion, a small amount of energy AE is in¬ 
jected at constant photon number into a blackbody of refer¬ 
ence temperature Tq. Once AE is fully comptonized, we are 
lef t with a Bose-Einstein spect rum with chemical potential 
llSunvaev fc Zeldovich|[l970cl l. The intensity difference is 



Here, x = hv/k^To, ks and h are the Boltzmann and 
Planck constants, and G{x) = xe^/{e^ — 1)^ describes a 
pure temperature shift. Note also that 7r^/[18()(3)] « 0.4561. 
The //-distortion is generated by inefficient diffusion of 
the photons in energy thro ugh scattering off of electrons 
JZeldovich fc Sunvaevlll963 1. The change in intensity reads 



We must furthermore take into acco unt the uncertainty 
in the temperature To = 2.726K (iFixsen et al.l 119961 : 

lFixsenll2009h of the refere nce blackbody up to second order 
dGhluba fc Sunvaevll20o3l . At this order, the deviation from 
the reference blackbody intensity i s the sum of a pure te m- 
perature shift and a //-distortions (IGhluba fc Jeon3l2014l ~). 


I^iu) = 


2hi/® 


G(a:)A(l + A) + y(a:)^' 


( 3 ) 


Hence, a relative error of A = 5 x 10”“^ = ATq/Tq in the 
precise value of To generates a //-distortion of amplitude 
y ~ 10“^. Patch-to-patch fluctuations in the GMB tem- 
perature will also contrib ute an average temperature shift 
dGhluba fc SunvaevI |20 o 3 1. which can simply be absorbed 
into the variable A. 

The recombination spectrum, /'"“(zz), can be eval¬ 
uated numerically with high accuracy. Since we are 
only interested in a detection of the recombination 
spectrum, we will use the template obtained from the 
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Detecting the recombination signal 3 


computation ol'^l Rubino-Martm. Chluba & SunvaevI 

([2OO6II: Ghluba & SunvaevI (2006ll: 

Ghluba. Rubino-Martm & Sunvaev 

(|200/tl and 

Rubino-Martm. Ghluba & Sunvaev 

(|200q). These cal- 


culations include the contributions from both hydrogen and 
helium. Refinements caus ed by helium feedback processes 
dChluba fc Sunva^ boid 'l are omitted here, but should 
not affect the main conclusions of this work. Similarly, 
small changes in recombination radiation at low frequencies 
(i/ < 1 GHz) caused by recombinations to highly exc ited 
states dChluba. Vasil fc Dursil[201(]l : I Ali-HaimoudI 1201^ 1 are 
omitted. 


dependence and adopt a value of 3x 10^ ^ Mp)/h broadly con - 
sistent with the best-fit HOD models of IZehavi et al.l d201lll . 
In all subsequent calculations, w e use the ha lo and sub- 
halo mass functi o ns pr ovided by [Tinker et al.l d2008fl and 
[Tinker fc Wetzel d2010[i . and integrate the sub-halo mass 
function from a minimum halo mass Mmin = Mq/Zi to 
the parent halo mass M. 

We assume that the same luminosity-mass relation 
holds for both central and satellite galaxies, and relate the 
galaxy infrared luminos ity to the host ha lo mass through 
the parametric relation dShang et al.[[2oi^ l 

i(i+.).(M, z) = Lo$(^)E(M)e[(l + z)u] , (9) 


2.2 Halo model of the GIB 

Our model of the GIB is based on a halo model description of 
the relation between star-forming galaxies and dark matter 


where Lq is an overall normalisation that must be con¬ 
strained from measurement of the GIB specific intensity 
in a given frequency range. The term 'I>(2) describes the 
redshift-dependence of the normalisation. We adopt the 
power-law scaling, = (1 -|- z)^'^, whereas, for the de- 


)f[Knox et al.l ([200 ill and 

Amblard & Gooravl ([200711 pendence S(M) of the galaxy luminosity on halo mass, 

olDe Bernardis & Gooravl 

201'^ for a model based on assume a log-normal distribution with mean mass 


The average GIB brightness at a given frequency (in 
unit of Jy sr~^) is 






(4) 


where x(D the line-of-sight comoving distance to redshift 
2 and 

L, 


j^{z) = dLn{L,z) 


dvr 


(5) 


is the mean emissivity of galaxies below a certain flux 
limit at a frequency v and per comoving volume . Here, 
and n{L, z) denote the infrared galaxy luminosity 
(in W Hz“^) and galaxy luminosity function, respectively, 
while (l + z)v desig nates the rest-frame frequency. Following 
[Shang et al.[ ([201 3), we split the mean emissivity into a sum 
of two contributions, 


in [Planck Gollaboration et al.[[2014lf) . Mefi characterises the 
peak and um/l the range of halo mass that produces a given 
luminosity L. Although our choice of = 0.5 is some¬ 

what arbitrary, the GIB angul ar power spectra t urn out to 
be fairly insensitive to um/l ([Shang et al.[[2M2 l. In addi¬ 
tion, the value of Meff is consistent with the peak of the 
stellar-to-halo mas s ratio inferred fr om semi-analytic galaxy 
formation models |[Guo et al.[[201(]h . 

For the galaxy spectral energy distribution (SED), we 
assume a modified blackbod y shape with a power-law emis¬ 
sivity as in [Hall et ^ ([2010[ i , 

u^B„{Td) V < Vo 

v~^ V > Vo ’ 


Q{v, z) oc 


( 10 ) 


where Bi,{T) is the brightness of a blackbody with tempera¬ 
ture T, and Td is the dust temperature. The grey-body and 
power-law connect smoothly at vo provided that 


Mz) = ldM^{M,z)[DiM,z) + DiM,z)] , ( 6 ) 

where 

D{M, z) = ^ArcL,,(i+,),,(M, z) (7) 

ftiM,z) = ^ Jdm^{M,z)L^^(i+^)^{m,z) (8) 

are the average emissivity produced by the central and satel¬ 
lite galaxies of a given halo at redshift z, dN/dM and dn/dm 
are the halo and sub-halo mass functions, and M and m 
are the parent halo and sub-halo masses. The numbers Nc 
and Ns of central and satellite galaxies are specified by the 
halo occupation distrib ution (HOD) (|Berlind fc Weinberg[ 
[2OO2I : [Zheng et al.[[2005ll . The total number of galaxies in 
a given halo of mass M thus is + Ng. Numerical simula- 
tions indicate that N c typically follows a step-like function 
([Kravtsov et al.[|2004 l. while Ns can be parametrised by a 
power-law with logarithmic slope a « 1. For the characteris¬ 
tic mass Mcen of the step function, we ignore any luminosity 


^ The data is available at: 

http://vivaldi.11.iac.es/galeria/jalberto/recomb/ 


dln0(i', z) 


d\nv 


= -7 


( 11 ) 


Unlike ISerra et ^ ([2014[ i. however, we will treat both the 
dust emissivity index /3 and the average dust temperature Td 
as free parameters, while we fix the SED index gamma to its 
mean v alue of 7 = 1.7 found by [Planck Gollaboration et all 
([2014b[l . The reason is that, at frequencies v < 10^ GHz, 
a change in the value of /3 tilts the GIB intensity in a 
way that could mimic a /r/y-distortion or the broadband 
shape of the recombination spectrum. This is illustrated 
in Fig[T] in which we show logarithmic derivatives of the 
average GIB intensity w.r.t. our free model parameters 
(Mefj, Td, ft, To) around the fiducial value (Mefj, Td,P,Lo) — 
(4.43 X 10^^ Mo/Zi, 26 K, 1.7, 4.8 x 10"''° W Hz'^). Glearly, 
the logarithmic derivative 91n7°'®/(?/3 cannot be expressed 
as linear combination of the other derivatives, so we must 
treat P as free parameter since it is not well constrained by 
the data. We will adopt a fiducial value of /? = 1.7. 

Note that the spectrum of star-forming galaxies typ¬ 
ically features broad emission lines in mid-infrared (from 
10 to 30/rm) caused mainly by po lycyclic aromat i c hy¬ 
drocarbon (PAH) molecules ( e-K- Lege r fc Pug ^ [1984 
[Allamandola. Tielens fc Barked [l985 : Smith et al. [2007 1. 
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4 Desjacques et al. 



Figure 1. Logarithmic derivatives /dpaii') of the average 

CIB intensity w.r.t. the model parameters p = Lq). 

The results are normalised such that all the derivatives are equal 
to unity for u = 0.1 GHz. Dotted curves indicate negative values. 


but these are far on the blue side of the spectrum (z/ > 
lO"* GHz). 


2.3 CIB anisotropies 

While we will ignore fluctuations in the and y-distortions 
across the skjQ, we must take into account anisotropies 
in the CMB and CIB intensity. In principle, it should be 
possible to subtract at least partly both CMB and CIB 
anisotropies from the observed patch using, e.g., data from 
the Planck experiment. In our fiducial analysis, however, 
we will assume they have not been removed. Therefore, we 
must include them in the Fisher information as they will 
contribute to the covariance of the signal, especially if the 
surveyed patch is small. 

The CIB intensity at a given frequency p and in a given 
direction fi can be expressed as the line of sight integral 



a{z)j^{z) 


A , Sj^{x{z)n,z) \ 

V Mz) J' 


Using the Limber approximation llLimbeilll954 l. the angu- 
lar cross-power spectrum of CIB anisotropies at observed 


^ The largest fluctuations are related to the spatially-varying y 
distortion signal caused by warm-hot intergalactic medium and 
unresolved clusters (e.g., I(5en Sz OstrikeJll99Sl: lRefregier & et aP 

I2QO0I : iMiniati et al.l I2OO0I : Izhang^Pe^^^^^ l 20 o 3 ) . while sig- 

nificant spatially-varying p distortions can only be created by 
anisotropic energy release processes. 


frequencies u and u' is (|Knox et al. IH) 

A) = ^ ^ a^{z)j^{z)j^,{z) (12) 

xP^^'^\k = £lx,z) . 

Assuming that spatial variations in the emissivity trace fluc¬ 
tuations in the galaxy number density, i.e. 5j/j = Sngjng, 
the power spectrum Pj is equal to the galaxy power spec¬ 
trum Pg. In the framework of the halo model discussed 
above, we thus have 

z) = PIA"(fc, z) + Pj^"(fc, z) + shot-noise . (13) 

Hereafter, we will ignore the shot-noise contribution because 
it is typically smaller than the 2-halo term on large scales. 
The contributions from galaxy pairs in the same halo, the 
so-called 1-halo term, is 

Pm ^ fdM ^\f:{M,z)f:,{M,zMk,z\M) 

jujv' J dM L 

+ f^r{M,z)fiiM,z)u{k,z\M) 

+ fl{M, z)fl, (M, z)u\k, z\M)\ , (14) 

where u(k,z\M) is the normalised Fourier transform of the 
halo density profile, whereas the 2-halo term reads 

P2^%,z) = ^D4k,z)D,,{k,z)Pun{k,z) , (15) 

JvJv' 

with 

D,{k,z) = jdM^h^{M,z)u{k,z\M) (16) 

X [rAM,z) + fl{M,z)\ . 

Here, Pun is the linear mass power spectrum extrapolated to 
redshift z, whereas 6i is the (Eulerian) linear halo bias that 
follows from a peak-background split applied to the halo 
mass function d N/dM. For consiste ncy, we use the fitting 
formula given in [Tinker et al.l 1I2OIOI I. The 1-halo term can 
only be seen with high angular resolution surveys reaching 
out to arcmin scales. 

The CIB intensity reported in each pixel of the surveyed 
patch can be written as 

/°™(j/,hi) = JdhI^^^{p,n)B{h,ni) (17) 

where B(n, n;) is the beam profile. Specialising to az- 
imuthally symmetric beam patterns, P(h, hi) is a function 
of h ■ hi solely, and can thus be expanded in the Legendre 
polynomials Pi'. 

1 °° 

P(h . hi) = — y {2£ + 1) BePiin ■ h,) . (18) 

dvr ' 

While this is an approximation that is not met in most cur¬ 
rent experimental designs, it is straightforward to generalise 
the computation. Our main conclusions will not change. For 
a single Gaussian beam profile, 

with h • hi = cosjI and ctb is related to the full width 
half maximum i9b of the beam in radians through rls « 
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2.3548ctb- I n this case, the series c oefficients are well approx¬ 
imated by (ISilk fc Wils^ 1 19801 : iBond fc Efstathiouf Il983 : 
IWhite fc Sredmckilll995h 

Bi « exp . (20) 


The CIB intensity averaged over the surveyed patch reads 
^ j dfii JdnD^^{v,n)B{h--hi) (21) 

= f dfii fdnD^^(u,h)YD*{n)YD{ni) ■ 

® /9™ fl ~ J 


Here, Hs = dvr/sky is the area of the surveyed patch. For an 
azimuthally symmetric patch centred at d = 0 (for simplifi¬ 
cation), the integral over the unit vector hi trivially is 



dhiYrih,) = WlSmO. 


( 22 ) 


For a circular cap (i.e., unit weight for all pixels with polar 
angle such that i9 < t) w), the harmonic transform of th e 
survey window is (e.g, iManzotti, Hu fc Benoit-Levvll20l3 l 


We 


y^[A+i(a:)-P,_r(a:)] . 


Here, Pe{x) denote Legendre polynomials and x = cosdw = 
1 — 2/sky is the cosine of the opening angle. Substituting 
the multipole expansions for the survey mask and the beam 
profile, the partial sky average of 7'^’® {u, fi) simplifies to 

7(!z) = ^ WeBe jdhD^\u, h)Y^{h) . (23) 


In the linear bias approximation, the Fourier 
modes of the CIB fluctuations are determined by 
5jv{k,z) = jv(z)b{k, u, z)S{k, z). If we now successively 
write D 8-® Ill® Fourier transform of Sj^{k,z), 

substitute the plane-wave expansion and use the Limber 
approximation, then the frequency-dependent multipoles, 
can eventu ally be written as the line of sight integral 
(|Knox et al.ll200l] l 



X a{z)D{z)Mz)^^!^^S{k)Ye°ik) . 


(24) 


Here, D{z) is the linear growth rate and z is defined through 
the relation kx{z) ~ i. The monopole of the CIB intensity 
computed from patches of the sky with coverage fraction 
/sky = Hs/(47r) is thus given by 

t{v)^^y^WeBeaei,{v), (25) 


where the multipole coefficients aiQ{i') are given by Eq. (|24l). 
Of course, the average (7‘^'®(j/)) over all such patches is 
exactly given by WoBoaoo{ie) = However, owing 

to cosmic variance, the monopole fluctuates from patch to 
patch with an amplitude given by 

[5D^^{v)Y = (26) 

= 772 I] WeWl,BiBl,{aeo{v)ae,o{u)) 







Figure 2. Contribution of galactic foregrounds and LSS-induced 
y-distortions to the intensity as a function of frequency. The var¬ 
ious curves represent the CMB monopole, galactic dust and syn¬ 
chrotron emission. We have also shown the monopole of the y- 
distortion (primordial -h LSS) with amplitude y = 5 x 10“^, as 
well as its expected rms variance (generated by LSS thermal SZ ef¬ 
fect) for a nearly all-sky survey (/sky = 0-8) with a beam FWHM 
of 1.6 deg. The specific intensity is in unit of Jansky per unit solid 
angle, Jy sr~^, where 1 Jy = 10“^® W m“^ Hz“^. 


where the cross-power spect rum lA _ is given by 

Ea. (|12[l . Note that, following iKnox et al.l (|2Q0lh . we have 
used the Limber approximation and approximated the ratio 
of Gamma functions squared in Eq. as « 2/7 to simplify 
the ensemble average (^ato{v)al,Q{v)). Therefore, our result 
strictly holds for 7 1. 


2.4 CMB primary and secondary anisotropies 

CMB anisotropies will also contribute to the covariance of 
the signal if they are not taken out. Ignoring the weak fre- 
quency dep endence induced by Rayleigh scattering (e.g., 
lLewisll2ni3h . the residual CMB intensity (i.e., with the ref¬ 
erence blackbody spectrum subtracted) averaged over the 
surveyed patch reads 

7°“®(i/) « 7^(i/) with TssA-te. (27) 

Here, A is the relative uncertainty on the temperature of the 
monopole and 0 is the average temperature anisotropy in 
the surveyed window. In analogy with the CIB, the variance 
of CMB intensity fluctuations across different patches of the 
sky reads 

[5P^^{v)f = ^Y.\WeBe\^Cf^\v,v) , (28) 

® l>2 

where the CMB intensity power spectrum is given by 

G\x)Ce. (29) 
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The contribution from Thomson scattering takes on the 
standard expression: 

2 

Ci = - dkk^gli{k)P^{k), (30) 

Jo 

where greik) denotes the pho ton transfer function, which 
can be obtained from camb (jLewis. Challinor fc Lasenbvl 
l2000li . We will assume that CMB intensity fluctuations are 
fully characterised by the primordial scalar amplitude As- 

The variance of primordial g- and ^/-distortions fluc¬ 
tuations (created before recombination) over patches can 
be ignored, since the current limits on the magnitude 
of y and g are fairly small so that fluctuations are ex¬ 
pected to contribute at the SI^/Iv — 10“^° level. Larger 
fluctuations of g and y across the sky could be cre¬ 
ated by the dissipation of acoustic modes with modu¬ 
lated small-sca le power due to non-Gaussianity in the ultra- 
squee zed limit llPaier fc Zaldarriagall2of^ : lGanc fc Komatsul 
l2012li. however, current upper lim its on /nl 2.5 ± 5.7 
llPlanck Collaboration et al.ll20153 l suggest that this case is 
unlikely for scale-invariant non-Gaussianity, so that we ig¬ 
nore it here. 

Secondary //-distortions arise in the late-time Universe 
because CMB photons scatter off hot e lectrons present in the 
gas of filaments and galaxy clusters llSunvaev fc Zeldovichl 
Il970bl i. The magnitude of this y distortion is an integral of 
the electron pressure along sight lines passing through the 
large scale structure (LSS). In a typical cluster, the electron 
temperature is Te ~ 1 keV and leads to i/ ~ 10“®. Exist¬ 
ing catalogues of galaxy clusters can be used to remove at 
least part of the signal generated at low redshift (i.e., large 
angular scales). Nevertheless, we will be left with a resid¬ 
ual monopole, which can be absorbed into P{v) as long 
as one is not interested in the primordial //-distortion, and 
a fluctuating contribution, whose angula r power spectrum 
is (within the Limber a pproximation, see lPersi et al.l[l995l : 
iR.efregier fc et al1l200fll l 



where the constant yo is ~ 1.7 x 10“^® Mpc“^ for 

a helium fraction and baryon density consistent with Big- 
Bang Nucleosynthesis constraints, Tp is the volume-average, 
density-weighted gas temperature and Pp(k,z) is the 3- 
dimensional power spectrum of pressure fluctuations. 

For simplicity, we will assume that Pp(k,z) = 
bpPiinik, z) is a biased version of the linear mass power spec¬ 
trum. We adopt bp = 85 , which yields a predictio n consistent 
with the simulations of lRefregier fc et al.l (l2000l i and, thus, 
provides a realistic upper limit to the signal. We compute 
Tp as the halo virial temperature Tvir weighted by the mass 
function. 


fp = — fdMM^{M,z)n,,{M,z) , (32) 

pm J dm 

where pm is the present-day average matter density, and 
rvu(M,2) =6.03x ^ ^ ^ ^ 


IQis Mo//i 


X (1 -I- «) flm K . 


(33) 


Hence, Tp is a strongly decreasing function of redshift. We 
adopt Avir = 200 at all redshift to be consistent with our 
choice of the halo mas s func tion, and assume /? = 2/3 as in 
iKomatsu fc Kitavamal lll999li . Furthermore, we only include 
halos in the mass range Mmin < M < Mmax(z). While we 
consider a fixed minimum mass Mmin = 10^® M 0//1 0,we 
allow Mina.x{z) to Vary with redshift to account for the pos¬ 
sibility of removing the SZ contribution from low redshift 
clusters. 

In Fig. [ 2 I we show the square-root of 

[SP^ i,y)] " = ( 34 ) 

“ e>2 

for an experiment with /sky = 0.8 and 1.6 deg angular 
resolution. The upper limit of the shaded region assumes 
Mmax = 10^® Mq//i at all redshifts. This corresponds to the 
signal that would be measured, had we not attempted to 
remove some of the //-distortion using external galaxy clus¬ 
ter catalogues. By contrast, the lower limit of the shaded 
region assumes that the SZ effect from groups and clusters 
with Tvir > 1 keV has been removed, so that Mmax[z) truly 
depends on redshift in our calculation. Note that these two 
limiting cases differ by < 50%, which shows that the contri¬ 
bution from virialized structures in filaments is quite sub¬ 
stantial. Overall, for the sky coverage adopted here, fluctu¬ 
ations in the //-distortions from the thermal SZ effect are of 
magnitude comparable to the recombination spectrum. Of 
course, [5/®^(:/)] is exactly zero for an all-sky survey. For a 
ground-based survey targeting a 16 deg^ patch of the sky as 
discussed below (see 113.21) . [5/®^(i/)] is approximately ten 
times larger than shown in Fig[2] 

We will hereafter ignore the signal covariance induced 
by the LSS through the thermal SZ effect as we will either 
consider satellite missions with /sky « 1, or a ground-based 
experiment for which (zz)] ^ and [(57‘^^®(zz)] ^ dominate 

the signal covariance. Notwithstanding, one should bear in 
mind that the LSS //-distortion could contribute significantly 
to the covariance as well as bias the primordial y interpre¬ 
tation if the sky fraction is significantly less than unity. 

For comparison, we also show in Fig[2] the average y- 
distortion for // = 5 x 10“^, together with the approximate 
level of galactic emission from synchrotron and thermal dust 
(for which we assumed a power-law and grey-body spec¬ 
trum, respectively). While we can reasonably handle the 
synchrotron emission with the low frequencies, removing the 
galactic dust monopole is more challenging as it is typi¬ 
cally brighter than the GIB, and has a similar grey-body 
spectrum. Furthermore, the galactic dust emission fluctu¬ 
ates significantly at large angular scales {I < 100), with a 
power spectrum Cg, oc (rather than Cg on i for the 

GIB). 


3 FORECAST FOR THE RECOMBINATION 
SPECTRUM 

Our model of the total (monopole) intensity Itot{v) mea¬ 
sured over a patch of the sky with coverage fraction f^ky 


® We have found that our predictions hardly change if we set 
Mmin = 10^^ Mq/Zj. 
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Detecting the recombination signal 7 


and in a frequency bin centred at i>a is 

+ D{vo.) 

+ r‘\v^) + D^\voi)+N{v^) , 


(35) 


where V°{ua) is the recombination signal which we wish 
to detect and the noise, N{va), includes contribution from 
the incident radiation and from the detector. For simplicity, 
we have assumed that foreground emission within our own 
Galaxy from synchrotron radiation of cosmic ray electrons 
and thermal emission from dust grains has been separated 
out from the signal using low- and high-frequency channels. 
We leave a more detailed analysis of foreground removal 
(including the CIB signal) for future work. 

The various signal components are shown in Fig. [S] to¬ 
gether with the f requency coverage and sensiti vity expected 
for the PIXIE llKoeut et al. 

|2014 1 and MILLIMETRON ( Smirnov et al.l 201^ sate llites, 
and a ground-based CII experiment ( Gong et al.ll2012i . The 
dotted-dashed curves represent the r.m.s. variance of GMB 
and CIB fluctuations in a patch of 16 deg^ assuming a (Gaus¬ 
sian) beam FWHM of 0.5 arcmin. For such a small sky frac¬ 
tion (/sky ~ 4 X 10“^), these are a few order of magnitudes 
larger than the signal in the frequency range 100-1000 GHz. 


3.1 Fisher matrix 


Consider Nv uncorrelated frequency bins centred at Ua and 
of bandwidth Sv covering the range (iZmin, r'max). Follow¬ 
ing l^gmarkj_IhyIor_&JHeaye^ (IIOOTII . the Fisher matrix 
is given bjQ 


Fij 



d{xc) i djxg) 
dpi “ dpj 


where (xa) = (l*°\va)) is the ensemble average of the resid¬ 
ual intensity over random realisations of the surveyed patch, 
Ca is the covariance of P°\ua), and p is the vector of model 
parameters. In our model, the signal ensemble average is 

(Xa) = P{Va) + + /”°(!Za) -f D{Va) + , (36) 

whereas the covariance of the data is given by 

Cs. = + (Af"(j/a)) (37) 

= [cf lyc) + !z.)] -b &ll . 

“ €>1 


For an all-sky survey, the covariance is simply Ca = SD. 
Furthermore, we ignore dependence on cosmological pa¬ 
rameters and assume that the recombination signal is per¬ 
fectly known, so that we can perform an idealised cross¬ 
correlation analysis. The model parameters therefore are 
p = (/i, y, A, As, Meff, Td, ft, Lo,Ar), where Ar is the ampli¬ 
tude of the recombination signal which we seek to constrain. 
Our fiducial parameter values are fi = 10“^, y = 5 x 10“^, 
A = 5 X 10“^ and As = 2.2 x 10“®, together with CIB pa¬ 
rameters Mefj = 4.43 X 10^^ M 0 // 1 , Td = 26 K, P = 1.7 and 
Lo = 4.8 X 10"'‘° W Hz-^ 


^ Although the trace is unnecessary here because the signal is 
one-dimensional, we keep it for the sake of generality. 



Figure 3. Average intensity as a function of frequency. The var¬ 
ious curves represent the CMB and CIB monopole, the T, /r and 
y-distortions and the recombination spectrum. The dotted-dashed 
curves labelled as y /and \/^-^cib indicate the r.m.s. vari¬ 
ance of CMB and CIB fluctuations measured in a 16 deg^ patch 
of the sky with a beam FWHM of 0.5 arcmin. The horizontal dot¬ 
ted lines show the frequency coverage and sensitivity expected for 
the PIXIE, PRISM and MILLIMETRON space experiments, as 
well as for a ground-based CII experiment. Note that the galactic 
dust monopole (not shown on this figure) is a grey-body similar 
to, though usually brighter than, the CIB (see Fig. Ell. 


3.2 Experimental setup 

The sensitivity of an ideal instrument is fundamentally lim¬ 
ited by the noise of incident photons. This is already the case 
of the best bolometric detectors. For photon noise-limited 
detectors, a gain in sensitivity can only be achieved by col¬ 
lecting more photons through an increase in the number of 
detectors, collecting area and/or integration time. 

For the moment, we assume M{v) does not exhibit fre¬ 
quency correlations and is normally distributed with a vari¬ 
ance (V^) = which depends only on the bandwidth 5v. 
The analysis is performed assuming top-hat frequency fil¬ 
ters, so that the filter response is 1 if £ [va — D, Ua -I- 
Su/2] and zero otherwise. Regarding the frequency coverage, 
it is difficult to sample frequencies less than ~ 30 GHz from 
space, for which the wavelength becomes larger than the 
typical size of the device (the collecting area) developed to 
measure them. However, t his frequency regime can be tar- 
geted from the ground Isee lSathvanaravana Rao et al.ll2015l . 
for a recent analysis). Finally, the spectral resolution will be 
= 15 GHz for the PIXIE experiment, but it could be as 
high as 1 GHz for a satellite mission like MILLIMETRON. 
For the sake of comparison, we will also consider a ground- 
based CII mapping instrument which achieve sub-Gigahertz 
resolution and targets 0(100) GHz frequencies correspond- 
ing to CII fine structure lines. We refer the reader to 
ISathvanaravana Rao et Tl (l2015h for ground-based surveys 
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8 Desjacques et al. 


Table 1. Experimental setups considered for the Fisher matrix forecast. The sensitivity of the ground-based experiment was derived for 
a 10m aperture telescope at frequency 238 GHz (see text for details). 




1?B 

SIi, [Jy sr 1] 

5u [GHz] 

i^-coverage [GHz] 

PIXIE 

full sky 

1.6 deg 

5 

15 

30 - 600 

PRISM 

full sky 

1.6 deg 

0.5 

15 

30 - 600 

MILLIMETRON 

full sky 

3 arcmin 

10“® 

1 

100 - 1000 

GROUND CII 

16 deg^ 

0.5 arcmin 

0.35 

0.4 

185 - 310 



Figure 4. Signal-to-noise ratio for the detection of the recombina¬ 
tion spectrum as a function of spectral resolution after marginali¬ 
sation over the remaining model parameters. We consider two dif¬ 
ferent frequency coverage, 30 < < 600 and 30 < < 3000 GHz; 

as well as two different scaling with the number of frequency 
channels: oc N^, ' (standard square-root law detectors) and 
oc (Fourier transform detectors). The SNR was normalised 

such that the two different types of detectors yield the same SNR 
for hv = 15 GHz. 


probing the < 1 GHz frequency range. Table [T] summarises 
the sensitivities of the various setups we consider in our 
Fisher forecast. Note that the lowest frequency measured 
is 30 GHz. 

We estimate the sensitivity of the ground-based exper¬ 
iment from the requirements given in iGong et al.l ll2012l) . 
Specifically, we consider a single-dish experiment with aper¬ 
ture diameter D = 10m. At 238 GHz (which corresponds to 
GII emission line at redshift z = 7 t hrough the transition 
^Pi /2 -f 2 Pi/ 2, see lGong et al.||20ll ), the resulting beam 
FWHM is approximately = 1.22A/D « 0.53 arcmin. The 
equivalent beam area is Gb ~ 8.82 x 10“® deg^. Assuming 
a noise equivalent flux density of NEFD=10 mjy s^^^ per 
spectral resolution element 5v = 0.4 GHz (consistent with 


the parameters in iGong et al.l [20121 ) . we can estimate the 
sensitivity from the relation 


NEFD 

■\/2TflB 


(38) 


where r is the total integration time and the factor of %/2 
arises from the Nyquist sampling and from the assumption of 
optical chopping. For a total integration time of 4000 hours, 
we find 5Ii, = 50 Jy sr“^ for a single detector. Averaging over 
20000 bolometers, we eventually obtain &Iv = 0.35 Jy sr“^. 

Let us assume that the frequency range is evenly split 
into frequency channels of spectral resolution Svq, and let SIq 
be the channel sensitivity that can be achieved at this reso¬ 
lution. Since the recombination spectrum features emission 
lines with a characteristic width of the order of their peak 
frequency (see Fig. [3)), we expect that the signal-to-noise 
ratio (SNR) for a detection of the recombination spectrum 
saturates below a certain spectral resolution. This, however, 
will be true only if the detector follows the standard square- 
root law SNR oc Ni, ' , where is the number of fre¬ 
quency channels corresponding to a spectral resolution Su. 
For a Fourier transform spectrometer (FTS), as utilised in 
PIXIE and PRISM, the s caling turns o ut to be SNR oc N~^. 

To see this, we follow IXogut et ^ (120111 ) and consider a 
total integration time r. If As is the number of time-ordered 
samples used for the Fourier transform, then each sample is 
observed during a time r/Na- Therefore, the noise in each 
time-ordered sample of the sky signal S{ti) is 


SS{ti) 


NEP 


(39) 


where the factor of two converts between time and frequency 
domains (the number of frequency channels is N^, = Na/2) 
and NEP is the noise equivalent power of the detected radi¬ 
ation. Note that NEP generally is a function of frequency. If 
the noise from different time measurement ti is uncorrelated, 
then in each frequency bin of the FTS it is 


5S{Ua,) 



2NEP 


(40) 


independent of the spectral resolution Su (at fixed Ni,5v). 
Therefore, since the sensitivity is 5Iv{va) oc 5S{va)/5v, we 
obtain the scaling 5Ii, oc Si'~^ oc N„. This implies that the 
signal-to-noise behaves like SNR oc 5I~^ oc N~^. This also 
means that (even under idealised conditions) averaging the 


® The source is measured one-half of the time, and one must 
differentiate the input signals. 
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Detecting the recombination signal 9 


SNR of individual synthesised FTS channels (improvement 
of ~ \/Nv) does not regain you the same sensitivity as di¬ 
rectly measuring at lower frequency resolution. 

In the FTS case, we thus expect the SNR to reach a 
maximum at some optimal spectral resolution before falling 
again as further decreases. Namely, under a change 
5vo —>■ in spectral resolution, the sensitivity of the de¬ 
tector becomes 


= Slo 



SIo 



(41) 


where N^Si, = NqSvo, while the signal-to-noise for, say, the 
model parameter pi scales like 


(SNR)^ 


1 

HI 


(1% 

\N^ 


2/3 N„ 

I E 



(42) 


Here, /3 = 1/2 and 1 for a square-root law and Fourier trans¬ 
form detector, respectively. Furthermore, we have assumed 
that both detectors yield the same SNR at the spectral res¬ 
olution Sl/Q. 

Figure |4] shows the SNR of the recombination spectrum 
as a function of spectral resolution for an all-sky satellite 
experiment that measures the intensity monopole, Eq. (1351) . 
with a sensitivity = 5 Jy sr“^ at spectral resolution 
5u = 15 GHz, very close to the PIXIE specifications. We 
consider two different frequency ranges extending up ei¬ 
ther to 600 GHz or 3 THz, so that the Lyman-a (Lya) 
line is only included in the latter case. The other compo¬ 
nents of our model [see Eq. (1351) 1 have all been marginalised 
over. Gonsider a standard square-root detector and mea¬ 
surements in the range 30 < < 600 GHz for instance. In 

this case, the SNR decreases mildly as the spectral resolu¬ 
tion increases from 10 to 120 GHz, and the Brackett-a (Bo?) 
line is gradually smoothed out. The SNR drops abruptly 
around 5v > 120 GHz, which corresponds to the disappear¬ 
ance of the Paschen-a (Pa) line. Therefore, one is basically 
left with only one line, the Balmer-a (Ha), at spectral reso¬ 
lution Sv > 120 GHz unless one extends the measurements 
out to ~ 3 THz to include the Lyman-a line (this is quite 
apparent in Fig. 0). In both cases, the SNR saturates for 
^ 10 GHz for a standard square-root detector whereas, 
for a FTS, the largest SNR is obtained for Sv ~ 50 and 
~ 200 GHz. Note that the spectral resolution of PIXIE is 
15 GHz, which is close to optimal for a standard square- 
root detector when the Lyman-a line is not measured. 


3.3 Mitigating the CIB contamination 

In order to assess the extent to which the CIB degrades 
the SNR of the recombination spectrum, consider a PIXIE 
experiment with sensitivity SI^, = 5 Jy sr~^, spectral resolu¬ 
tion Su = 15 GHz and frequency range 30 < < 600 GHz. 

Assuming we have perfect knowledge of all the model pa¬ 
rameters except the amplitude of the recombination signal, 
the SNR for the recombination spectrum is SNR ~ 0.44. On 
marginalising over the four parameters (Mes, Ti,/3, Lq) de¬ 
scribing the CIB, the SNR drops to 0.17. Further marginal¬ 
isation over the primordial CMB spectral distortions y and 
p brings the SNR down to 0.10. For the wider frequency 
coverage 30 < v < 3000 GHz, the degradation in the SNR 
following marginalisation over the CIB is nearly a factor of 



u [GHz] 


Figure 5. Effect of spectral resolution on the recombination spec¬ 
trum. At spectral resolution &u > 100 GHz, the Paschen-o line 
is wiped out and, unless the frequency coverage is wide enough 
to include the Lyman-a line, the SNR drops sharply (see Fig. 
0. The shaded region indicates the frequency range probed by 
the ground CII experiment. In this spectral region, the narrow 
feature at u ^ 270 GHz is the combination of two lines with 
positive and negati ve peak intensity imprinted during Hell— »HeI 
recombination Isee lRubino-Martm. Ghluba fc Sunvae\il200Sl l. 


4 (down from 0.96 to 0.26). Unsurprisingly, the amplitude 
of the recombination spectrum Ar strongly correlates with 
Meff, /3 and Lq (with a correlation coefficient jrj « 0.9 in all 
cases). This is another justification for treating /3 as free pa¬ 
rameter. This also highlights that the continuum part of the 
cosmological recombination radiation is more difficult to iso¬ 
late, even if its amplitude may be close to the sensitivity. Ac¬ 
cessing the variable component with its quasi-periodic fea¬ 
tures strongly improves the situation. For an experiment like 
PRISM, the recombination spectrum could be detected with 
a signal-to-noise ratio of SNR~ 2.4 if the frequency cov¬ 
erage extends up to 3 THz. Finally, for the MILLIMETRON 
experiment, the SNR is as large as ~ 2500 for the parameter 
values quoted in Table [H 

To ascertain the extent to which our CIB model pa¬ 
rameters could be constrained by a measurement of the 
CIB anisotropies , we consider the Fisher matrix (e.g., 
iPenin et aI.ll2C)l 


/r.ciB 


fsky 

2 


^(2^+l)tr 

£ 



dc, 

dpi ^ dpj 


(43) 


Here, Ct is the full, Ni, x Ni, covariance matrix at a given 
multipole £. The multiplicative factor of /sky reflects the fact 
that, for partial sky coverage, modes with a given multipole 
are partially correlated, hence the variance of each measured 


© 0000 RAS, MNRAS OOO.ITIfT^ 
















10 Desjacques et al. 


Ce is larger. The matrix takes the form 

/ + Afp 

Cr= 

V 

where, for shorthand convenience, {ui,Vj) and 

A /7 = Afi{u) is the angular power spectrum of the instru¬ 
mental noise. Assuming uncorrelated and isotropic pixel 
noise, the la tter is given by 

l|Knoxl flQQRh . where fipix is the solid angle subtended by 
one pixel, Upi* is the noise per pixel and Bi is the exper¬ 
imental beam profile, Eq. (HSl)- Since our convention is to 
express the brightness I(i') in unit of Jy sr“^, the CIB an¬ 
gular power spectra are in unit of Jy^ sr“^. The noise 

per pixel squared is Npi^SI^ whereas the pixel solid angle is 
Aix/Npi^. Therefore, the weight is simply = AirSIp 
and, as expected, is independent of the sky pixelization. For 
the PIXIE- and PRISM-like specifications, we find = 314 
and 3.14 Jy^ sr~^, respectively. For illustration, Fig.[6]shows 
the CIB angular power spectrum at two different frequencies 
for a PRISM-like experiment with 3 arcmin angular resolu¬ 
tion of the imager. 

We have evaluated the Fisher matrix Eq. (HHl) for these 
experimental setups with a conservative value of /sky = 0.8, 
and subsequently computed the SNR for the recombina¬ 
tion spectrum upon adding the Fisher information, Fij = 
-I- where is the Fisher matrix for a mea¬ 

surement of the intensity monopole. For our fiducial beam 
FWHM of 1.6 deg (with /max = 1000), the improvement 
is negligible. A measurement of the CIB anisotropies at a 
much finer angular resolution of 3 arcmin (for which we 
adopt /max = 10“*) somewhat increases the SNR, but the im¬ 
provement remains small. Namely, for the PIXIE and PRISM- 
like experiments respectively, the SNR for a measurement of 
the recombination spectrum increases from 0.10 to 0.16 and 
from 1.02 to 1.15 when the angular resolution is changed 
from 1.6 deg to 3 arcmin. This highlights that, in terms of 
the recombination signal, increased angular resolution is not 
a main issue. We have not pushed the angular resolution 
below 3 arcmin as our model does not include shot noise, 
which we expect to contribute significantly to the C^'' at 
multipoles £ > 10*. 

We have not computed F)'/*’ for the small field CII 
survey because the numerical evaluation is very time- 
consuming, owing to the exquisite frequency resolution. Nev¬ 
ertheless, Fig.[3]clearly shows that, since patch-to-patch fluc¬ 
tuations in the CIB are a few orders of magnitude larger than 
the recombination signal at frequencies ~ 200 GHz, remov¬ 
ing the CIB is an even greater challenge for ground-based 
surveys targeting those frequencies. In the idealised situa¬ 
tion where the CIB is perfectly known and the patch has 
been cleaned from CMB anisotropies, the signal-to-noise ra¬ 
tio is SNR ~ 2.74. The feature at ~ 270 GHz (see Fig. [S]) 
helps distinguish the recombination spectrum from the CMB 
spectral distortions. Overall, our Fisher analysis shows that 
an all-sky survey covering a large frequency range at mod¬ 
erate frequency resolution Av ~ 10 GHz will perform much 
better than an experiment targeting a small patch of the sky 
at high spectral resolution {Au ~ 0.1 GHz). 



multipole 

Figure 6. CIB angular power spectrum as measured by an ex¬ 
periment with PRISM sensitivity and a beam angular resolution 
of 3 arcmin. For the spectral resolution considered, 38 X 38 cross¬ 
power spectra C^'' can be computed at each multipole. Notice 
the flattening at £ > 1000 caused by the 1-halo term. 


4 CONCLUSIONS 

We have presented a Fisher matrix forecast for the detection 
of the recombination line spectrum with future space exper¬ 
iments. The main caveat of our analysis is the assumption 
that galactic synchrotron and dust emission have been sep¬ 
arated out using low- and high-frequency channels. We have 
also assumed that most of the thermal SZ effect generated by 
the large-scale structure can be removed using external cat¬ 
alogues of galaxy clusters when it contributes significantly 
to the signal covariance (when /sky 1). Aside from that, 
our main findings can be summarised as follows: 

• The cosmic infrared background is the main contam¬ 
inant to the recombination signal. Detecting the recombi¬ 
nation lines requires sub-percent measurement of the CIB 
spectral shape at frequencies f 1 THz. The T, y and 
/r-distortions only weakly correlate with the recombination 
spectrum if the frequency coverage is large enough. 

• Adding information from the CIB angular power spec¬ 
trum does not greatly improve the SNR, even at arcmin an¬ 
gular resolution. Note that we have not considered the possi¬ 
bility of cross-correlating the CIB with large-scale structure 
to better constrain the CIB model parameters. 

• While a ground-based Cll-like experiment targeting a 
small patch of the sky cannot beat the CIB fluctuations, 
even with very high spectral resolution and exquisite sensi¬ 
tivity, an all-sky satellite mission can do this. 

• For an all-sky measurement in the frequency range 30 < 
v < 600 GHz, a spectral resolution Sv < 10 GHz is optimal 
at fixed SvSIp. For a Fourier transform spectrometer, the 
optimal spectral resolution is larger {5v ~ 50 GHz) and 



© 0000 RAS, MNRAS 000. [HO 










Detecting the recombination signal 11 


dependent on the frequency coverage (see Fig. II, owing to 
the noise-scaling SIi, oc 5u~^. 

• A future all-sky satellite mission with sensitivity SIv = 
0.1 Jy sr“^ and spectral resolution Su = 15 GHz can de¬ 
tect the recombination lines at 5a for frequency coverage 
30 < v < 600 GHz. If higher frequency channels are in¬ 
cluded (30 < V < 3000 GHz), a sensitivity of 51 v — 0.25 
Jy sr“^ is needed for a 5cr detection. An experiment with 
milli-Jansky sensitivity in frequency channels 5v = 1 GHz 
like MILLIMETRON may measure the recombination lines 
with a SNR of ~ 2500. 


A few more comments are in order. Firstly, the GIB 
fluctuations are produced by fluctuations in the distribu¬ 
tion of high-redshift galaxies and, therefore, should correlate 
strongly with the fluctuations measured in large scale struc¬ 
ture (LSS) surveys provided the latter are deep enough to 
resolve stellar masses of order M* > 10^°“^^ Mq/Zi. There¬ 
fore, it may be possible to remove the GIB substantially if 
we overlap with a deep LSS survey. 

Secondly, a major issue for the detection of the recombi¬ 
nation signals, and any of the primordial distortions really, 
will be the calibration. To separate different frequency de¬ 
pendent components, a calibration down to the level of the 
sensitivity is required. In this way, different channels can 
be compared and the frequency-dependent signals can be 
separated. For the recombination signal, it may be enough 
to achieve sufficient channel cross-calibration, since in con¬ 
trast to the primordial and y-distortion the signal is quite 
variable. This issue will have to be addressed in the future. 


Finally, we highlight that a detection of the recombi¬ 
nation signal also guarantees a detection of the low-redshift 
^-distortion and the small-scale dissipation signals. These 
two signals are expected in the standard cosmological 
model and allow us to address interesting questions about 
the reionization and structure formation process, as well 
as the early Universe and inflation physics. Even with 
a sensitive low resolution CMB spectrometer it may 
furthermore be possible to transfer some of the absolute 
calibration to an independent high-resolution GMB im¬ 
ager. This could further open a possibility to extract the 
line-scatter ing signals from the dark ages and the recombi¬ 


nation era (BasUj_JIerrjandez 2 MQateMlJdflJL.SunYaey| 


Rubino-Martfn. Hernandez-Monteagudo fc SunvaevI 


2004 


20051: 


LewisI 2013l i. which in terms of sensitivity are also within 


reach. 


In summary, detection of the hydrogen and helium lines 
from recombination of the early universe is a hugely chal¬ 
lenging but highly rewarding goal for the future of CMB 
astronomy. This field, largely neglected for three decades, 
is ripe for exploitation. An all-sky experiment is required to 
measure the recombination lines. This can best be done from 
space, or possibly via long duration balloon flights. Likely 
rewards would include the first spectroscopic study of the 
very early universe and the first measurement of the pri¬ 
mordial helium abundance, as well as unsurpassed probes of 
new physics and astrophysics in a an entirely new window 
on the earliest epochs that we can ever directly access by 
astronomical probes. 
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